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ABSTRACT 


This technical memorandum documents the longwave radiation parameterization developed 
at the Climate and Radiation Branch, NASA/Goddard Space Flight Center, for a wide variety of 
weather and climate applications. Based on the 1996- version of the Air Force Geophysical 
Laboratory HITRAN data base (Rothman et al„ 1998), the parameterization includes the 
absorption due to major gaseous absorption (water vapor, C0 2 , O 3 ) and most of the minor trace 
gases (N 2 0, CH 4 . CFC's), as well as clouds and aerosols. The thermal infrared spectrum is 
divided into nine bands. To achieve a high degree of accuracy and speed, various approaches of 
computing the transmission function are applied to different spectral bands and gases. The 
gaseous transmission function is computed either using the ^-distribution method or the table 
look-up method. To include the effect of scattering due to clouds and aerosols, the optical 
thickness is scaled by the single-scattering albedo and asymmetry factor. The optical thickness, 
the single-scattering albedo, and the asymmetry factor of clouds are parameterized as functions of 
the ice and water content and the particle size. The parameterization can accurately compute 
fluxes to within 1% of the high spectral-resolution line-by-line calculations. The cooling rate can 
be accurately computed in the region extending from the surface to the 0.01 -hPa level. 

The computer code for this longwave parameterization is very easy to use. It has been 
implemented in atmospheric models for cloud, weather, and climate studies at many government 
institutes and universities. The code and sample calculations are accessible at 
http://climate. gsfc. nasa. gov/~chou/cli rad_lw. 
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1. INTRODUCTION 


Thermal infrared (or longwave, LW) radiation plays a crucial role in affecting weather, 
climate, and the sensitivity of climate to external radiative forcing. It is, therefore, important to 
have an accurate LW radiation parameterization in atmospheric models for weather and climate 
studies. In numerical model simulations of weather and climate, calculations of LW fluxes can 
easily take a third or more of the total computing time. As the spatial and temporal resolution of 
the models increases and the treatment of physical processes improves, it becomes clear that we 
need a fast and accurate LW radiation parameterization. 

Detailed calculation of the LW fluxes involves three sets of integration: spectral integration, 
vertical integration, and directional integration. The spectral integration is the most time 
consuming, but it is the vertical integration that makes the parameterization for the LW radiation 
particularly difficult. If the atmosphere were vertically homogeneous in pressure and 
temperature, then the absorption coefficient would be a function only of wavenumber, and 
wavenumbers with the same absorption coefficient would be radiatively identical. The spectral 
integration could then be greatly simplified by using the ^-distribution method (c.f Arking and 
Grossman, 1972), in which only one set of radiative transfer calculations is needed for all the 
wavenumbers with the same absorption coefficient. In the real atmosphere, however, the 
dependence of the absorption coefficient on pressure and temperature varies with wavenumber, 
and in principle no two spectral intervals can be treated identically in the radiative transfer 
calculations. The difficulty is augmented by the narrowness of molecular absorption lines, which 
makes the absorption coefficient vary rapidly with wavenumber. Thus, it requires a very high 
spectral resolution to obtain spectrally integrated fluxes with high accuracy. To further 
complicate the flux calculations, there are numerous absorbers that have to be taken into account. 
Except for water vapor, those absorbers spread over the LW spectrum in narrow bands. The 
effects of those bands need to be computed individually. 

At the Goddard Climate and Radiation Branch, we have developed various LW radiation 
parameterizations for the major water vapor, C0 2 and 0 3 absorption bands (Chou and Kouvaris, 
1991; Chou et al., 1993). These IR radiation parameterizations have been shown to be accurate 
and efficient in computing cooling rate not only in the troposphere and lower stratosphere, but 
also in the middle atmosphere (up to the 0.01 -hPa level). The code was documented in a NASA 
Technical Memorandum (Chou and Suarez, 1994). 
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This technical memorandum is an extension of Chou and Suarez (1994) by including the 
absorption due to minor trace gases of N 2 0 , CH 4 , and CFCs (Kratz et al., 1998) and 
parameterizations for the absorption and scattering due to clouds (Chou et ah, 1999). Aerosols are 
also included, but with optical properties specified as input parameters. This LW radiation code 
has been implemented in the Goddard general circulation models (Bacmeister and Suarez,, 2001; 
Sud and Mocko, 1999) and the cloud ensemble model (Tao et al., 2001). It has also been used at 
various universities and government institutes. 


2 . INFRARED TRANSFER EQUATIONS 

Let us consider a thin atmospheric layer at pressure p ' with a temperature 0', which has a 
differential pressure thickness dp' and an optical thickness dx v at wavenumber v (Figure 1). The 
radiance emitted by this layer is R v (W)dx w /p, where p is the cosine of the angle between the beam 
and the vertical, and R v is the Planck function. Let us further consider a lower level at pressure p 

at which we wish to compute the downwelling flux, and let us assume that between p and p' 
there is a mixture of gaseous, cloud, and aerosol absorbers with monochromatic optical thickness 
T'viPtP')’ The transmittance between p and p in the direction p is e~ x ^ p ' p ^ When scattering 

is neglected, the contribution of the layer between p' and p' +dp' to the radiance at the level p 
is then given by 

^/ v i (n) = /^(e')^ L e _Tv(/ '’ p ' ) ^ (2.i) 



P’ 

p' + dp' 


P 


Figure 1 . Contribution from a differential layer at p' to the downward radiance at p. 
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The contribution of the layer to the downward flux at p can be derived by angular 
integration of (2. 1 ) over the hemisphere, 

dF v l = 2n f 

Jo 

= 2nR,(Q') f (2.2) 

Jo P 

= -By (6 ) dT v (p,p') 

where B l/ (Q)=nR,(Q) is the Planck flux, and T v (p,p') is the flux transmittance for isotropic 
radiation given by 


7v(A/>') = 2 f c 

Jo 

The total downward flux integrated over heights and wavenumbers is then given by 



( 2 . 3 ) 


( 2 . 4 ) 


Similarly, the upward flux can be expressed as 

|e v ^ < e «)+ 0 - £ « K <(’■(>,) 


F (p)= \dv 




( 2 . 5 ) 


where ^ is the surface emissivity, p s is the surface pressure, and 0* is the Earth's surface skin 
temperature. Finally, the cooling rate is given by 


30 (p) 
3 1 



(2.6) 


where c p is the heat capacity of air at constant pressure, and g is the gravitational acceleration. 
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3. SPECTRAL BANDS 


3.1 The 10 Bands 


For computing thermal IR fluxes, the spectrum is divided into 9 bands and a subband. 
Table 1 shows the spectral ranges for these 10 bands, together with the absorbers involved in each 
band. The water vapor line absorption covers the entire IR spectrum, while the water vapor 

continuum absorption is included in the 540-1380 cm 1 spectral region. The absorption due to C0 2 
is included in the 540-800 cm 1 region, and the absorption due to 0 3 is included in the 980-1 100 
cm region. The minor absorption due to CH 4 , N 2 0, CFC’s , and C0 2 is scattered between 800 

cm 1 and 1380 cm ‘ region in Bands 4-7. The absorption due to N 2 0 in the \l-jjm region is 
included in Sub-band 3a and is identified as Band 10. 

3.2 Planck-Weighted Band Integrals 

From (2.4) and (2.5), the downward and upward fluxes integrated over a spectral band i can 
be written as 


F,hp) = 



(dT(p,p )) , 

\-*rr 


F? (p)=e i B i (d s )T i (p,p s )+ 



+ {\-e l )F i l (p s )T i (p,p s ) 


(3.1) 


(3.2) 


where 


F,(P)= 1 

» 

F v (p)d\ 

(3.3) 

J 

Av, 


B/0) = 

• 

Rj(Q)dv 

(3.4) 

• 

'av, 


Ti(p,p') 

f W K(p.p’)dv 

(3.5) 
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Table 1 . Spectral bands, absorbers, and transmittance parameterizations. 


Band 

Spectral Range 

(crrr i) 

Absorber 

Options for 

Transmittance 

Parameterization 

"LOW” "HIGH" 

1 

0-340 

H2O line 

K 

T 

2 

340-540 

H2O line 

K 

T 

3a 

540-620 

\ H2O line 

K 

T 

3b 

620-720 

f H2O continuum 

C 

C 

3c 

720-800 „ 

1 co 2 

K 

T 

4 

800-980 

H2O line 

K 

K 



H2O continuum 

C 

C 



C 0 2 , FI 1 , FI 2 , F 22 

K 

K 

5 

980-1 100 

H2O line , CO2, FI 1 

K 

K 



H2O continuum 

C 

c 



0 3 

T 

T 

6 

1100-1215 

H2O line 

K 

K 



H2O continuum 

C 

C 



H 2 0 , CH 4} F 12 , F 22 

K 

K 

7 

1215-1380 

H2O line , N 2 0 , CH 4 

K 

K 



H2O continuum 

C 

C 

8 

1380-1900 

H2O line 

K 

T 

9 

1900-3000 

H2O line 

K 

K 

10 

540-620 

H2O line , CO2 , N 2 0 

K 

K 



H2O continuum 

C 

C 


K : /c-distribution method with linear pressure scaling. 

T: Table look-up with temperature and pressure scaling. 

C: One-parameter temperature scaling for water vapor continuum absorption. 
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(3.6) 


dJXpjj ) 

dp' 



3r v (/7,p ) 
d~P " 


Jv 


are the Planck-weighted band integrals, and Av, is the width of the i th band. 


Within each spectral band, either the range of B v is sufficiently small or the shape of is 
sufficiently independent of temperature that we can make the following approximation 


T ‘(P'P' , f Bv(Q 0 )Tv(p,p' )dv 

5,-(0„) J Av, 


dT(p,p) 
dp' , 


B,(%) 


f Wo) 

Jav. 


dT v (p,p ' ) _ dTjjpj ^ ) 

d P ' d/7’ 


(3.7) 


(3.8) 


where 0 f> is a typical value of the atmospheric temperatures, which is set to be 250 K. 


With the approximation (3.7) and (3.8), the band-integrated fluxes reduce to 

F, r 1 (P ) = f 5,(0 ' ) dTl( £ P ' - dp' (3.9) 

Jo op 

F?{p) =e i B i (Q s )T i (p,p s )+ f%(0’) ^£Ml dp 

J Ps d P (3.10) 

+ ( l ~£ l )F' i l (Ps)Ti(p,p s ) 

3.3 Band- Integrated Planck Functions 

The spectrally integrated Planck fluxes were pre-computed for each band and then fit by a 
^-degree polynomial in temperature, 

5 

m) = Q.0 "b (3. 11) 

n = 1 

When integrated over all bands, errors in this regression are negligible (< 0.1%) for 150 K < 0 < 
350 K. The coefficients c’s are listed in Table 2. 
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Table 2. Coefficients for computing the spectrally integrated Planck flux from Equation 
(3.1 1). Units of 0 in (3.1 1) are Kelvin. 


Band 

c o 


C 2 

Cj 

c 4 

C S 

i 

5.344e+0 

-2.062e-l 

2.533e-3 

-6.863e-6 

1.012e-8 

-6.267e-12 

2 

2.715e+l 

-5.404e-l 

2.950e-3 

2.723e-7 

-9.338e-9 

9.968e-12 

3 

-3.486e+l 

1.11 3e+0 

- 1 .301 e-2 

6.496e-5 

-1.1 82e-7 

8.042e-l 1 

4 

-6.05 le+l 

1.409e+0 

-1.208e-2 

4.405e-5 

-5.674e-8 

2.566e-l 1 

5 

-2.669e+l 

5.283e-l 

-3.445e-3 

6.072e-6 

1.252e-8 

-2. 1 55e- 1 1 

6 

-6.727e+0 

4.226e-2 

1.044e-3 

-1.292e-5 

4.740e-8 

-4.486e- 1 1 

7 

1.879e+l 

-5.836e-l 

6.968e-3 

-3.939e-5 

1.012e-7 

-8.230e-l 1 

8 

1 ,034e+2 

-2.513e+0 

2.375e-2 

-1.069e-4 

2.184e-7 

-1.370e-10 

9 

-1.048e+l 

3.82 le-1 

-5.227e-3 

3.441e-5 

-1.108e-7 

1.409e-10 

10 

1.677e+0 

6.540e-2 

-1.81 3e-3 

l.291e-5 

-2.672e-8 

1.979e-l 1 


3.4 Sensitivity of Upward Flux to Surface Temperature 


The parameterization also computes the sensitivity of upward fluxes at all levels to the 
surface temperature 0 S . In a general circulation model, the LW parameterization is called 

relatively infrequently (every few hours), whereas the boundary layer and land surface 
parameterizations use time steps of a few minutes. To maintain consistency between the upward 
radiative fluxes aloft and that at the surface when the surface temperature changes significantly 
between calls to the LW parameterization, all fluxes are linearized about the surface temperature 
at the beginning of the radiation interval, and radiative heating rates are recomputed based on this 
linearization every few minutes, at each time step. The partial derivative of the upward flux with 
respect to the surface temperature is computed from 


dFjhp) 

30 , 


= Ti(p,p s ) 


e,3ff(0) 

30 


-te=e. 


Ti(p,p s )£i 



n= 2 


n - 1 

5 


(3.12) 
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4. BAND- AVERAGED TRANSMISSION FUNCTIONS OF GASES 


For accuracy and speed considerations, the flux transmittance defined by (3.7) is computed 
using three different approaches, depending on the absorber and the spectral band: 

• The ^-distribution method with linear pressure scaling is applied to the absorption due to 
water vapor and minor trace gases. This method is also applied to the C0 2 absorption in 
Band 3 if accurate cooling rate calculations in the middle atmosphere (p < 10 hPa) are not 
required. 

• The transmittances due to C0 2 and 0 3 absorption in Bands 3 and 5 are obtained from pre- 
computed transmittance tables. Because the /c-distribution method with linear pressure scaling 
underestimates the water vapor cooling rate in the middle atmosphere, the transmittances of 
the three strongest water vapor absorption bands (Bands 1, 2 and 8) are also derived from pre- 
computed transmittance tables if accurate computations of the water vapor cooling in the 
middle atmosphere are required. 

• The transmittances due to water vapor continuum absorption in Bands 3-7 are computed 
using a one-parameter scaling approach. 

Applications of these parameterizations to the different spectral bands and absorbers are 
summarized in Table 1 . 

4. 1 The k-distribution method 

As shown in Chou et al. (1993), the cooling due to water vapor in the lower atmosphere {p 
> 10 hPa) is primarily attributable to the spectral regions away from the center of absorption 
lines. The absorption coefficient is approximately linear in pressure, and its dependence on 
temperature varies smoothly with wavenumber. Under these conditions, the absorption 
coefficient at any temperature and pressure is simply proportional to its value at a reference 
pressure and temperature. The absorption coefficient can, therefore, be extrapolated from the 
absorption coefficient at a reference pressure and temperature according to 

K(p, Q) = K(p r , e j (p /p r ) m h re, e r ) (4. t ) 

where p r and 0, are the reference pressure and temperature, m is an empirical constant, and h(Q. 
0J is the temperature scaling factor, which we fit by 
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h(Q, ej=y+<x(e-e r )+|3(e- ef 


(4.2) 


The regression coefficients a and (5 are derived for each spectral band to minimize errors in flux 
transmittances. Details of the derivation of the coefficients a and P are given in Chou et al. 
(1993). Values of p n 0 n m, a, and P are given in Table 3 for scaling the water vapor and C0 2 
absorption coefficient. 


Table 3. Spectral bands, reference pressure and temperature, and the coefficients for the 
pressure and temperature scaling used in the ^-distribution method. Band 3 is divided into three 
sub-bands. 


Band 

h 2 o p r 

= 500 hPa 

0 r = 250 K 

co 2 


e r = 

250 K 

m 

ol(K') 

P (K 2 ) 

p r ( hPa ) 

m 

o.(K‘) 

P(^ 

1 

1.0 

.0021 

-1. 01e-5 





2 

1.0 

.0140 

5.57e-5 





3a 

1.0 

.0167 

8.54e-5 

300 

0.50 

.0179 

1.02e- 

3b 

1.0 

.0167 

8.54e-5 

30 

0.85 

.0042 

2.00e-: 

3c 

1.0 

.0167 

8.54e-5 

300 

0.50 

.0184 

1 . 1 2e- 

4 

1.0 

.0302 

2.96e-4 





5 

1.0 

.0307 

2.86e-4 





6 

0.77 

.0195 

l.lle-4 





7 

0.5 

.0152 

7.61e-5 





8 

1.0 

.0008 

-3.52e-6 





9 

1.0 

.0096 

1.64e-5 





10 

1.0 

.0149 

6.20e-5 

300 

0.50 

.0179 

1 .02e- 
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With the scaling of (4.1), the monochromatic flux transmittance (2.3) can be expressed as 
T v (p,p')=2\ ] (4.3) 

Jo 

where w is the scaled absorber amount given by 

1 CP 

w(p,p)=~ c{p")<p"p r ) m h{W'$ r )dp" (4.4) 

S J p* 

and c is the absorber mass concentration. It can be seen from (4.3) and (4.4) that the dependence 
of T v on wavenumber is separated from its dependence on pressure and temperature. 

Wavenumbers which have the same absorption coefficient at p r and 0 r will have the same 

transmittance at all other pressures and temperatures. Within a narrow spectral interval where the 
Planck function can be considered constant, these wavenumbers are radiatively identical, and the 
integration over wavenumber can be replaced by an integration over the absorption coefficient k , 

I ( )dv / Av = j ( )f(k)dk (4.5) 

Jav Jo 

where f(k) is the ^-distribution function at the reference pressure p r and temperature 0 r , and 

f(k)dk = 1 (4.6) 

This is the basic assumption of the ^-distribution method. It should be noticed that the 
approximation for relating the absorption coefficient at a pressure level to that at other levels, as 
shown in (4.1), is different from the correlated-k distribution method (e.g. Wang and Shi, 1988; 
Goody et al., 1989; Lacis and Oinas, 1991; Fu and Liou, 1992; Mlawer et al., 1997). However, it 
is also a type of the correlated-^ approximation with a well-defined correlation. 

If we divide a wide spectral band into a number of narrow sub-bands, for which the Planck 
function can be considered constant, and apply the ^-distribution method to each sub-band, the 
Planck-weighted transmission function given by (3.7) becomes 
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I 

Bj(Q 0 ) f T v (w)dv 

Jc,, 


j 

* , OV j 



B(e 0 ) 


I 

j 1 

Bj (6 0 )8v j 1 T k (w)fj(k)dk 

J k 

B(Q 0 ) 


where 


T k {w)=l \ e ( 4 - 8 ) 

Jo 

fl/ej is the mean Planck function of the sub-band j. and f/k) is the ^-distribution function of the 
sub-band at the reference pressure p r and temperature 0 r . 

The flux transmittance given by (4.8) can be accurately computed by using the following 
approximation (see Figure 5b of Chou and Arking, 1980) 

T k (w)=e- kwl » ( 4 . 9 ) 

where 1 /(I is the diffusivity factor taken to be 1.66. Using this approximation and replacing the 
k integral by a simple quadrate, Equation (4.7) reduces to 

N 

T( W ) = Ye~ k ^A gn (4-10) 

n= 1 

where A g„ is the Planck-weighted /(-distribution function for the m" 1 ^-interval given by 

= (4U) 

j 

Values of A g n are first derived from line-by-line calculations and then slightly adjusted using 

regression so that the rms difference between the transmittances computed from the line-by-line 
method and from (4.10) is minimized. (The line-by-line calculations are briefly discussed in 
Section Section 9). The adjustment is necessary because of the small number of terms (N = 6) 
used and the diffusivity (4.9) approximation applied. Details are given in Chou et al. (1993). 
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Calculations of the Planck-weighted flux transmission function using (4.10) are very fast 
for two reasons: (1) For each k , it requires only L exponential operations for calculating 

e / =1, 2, L, where L is the total number of atmospheric layers, and A w } is the scaled 

absorber amount in the layer /. The transmittance of a path consisting of more than one layer can 
be obtained by multiplying the transmittances of individual layers. By comparison, the 
commonly used band models require at least 0.5 L 2 exponential operations. (2) We choose values 
of k in such a way that 

K = l\Ki n=2,3, N (4.12) 

where rj is a positive integer. With this choice of k' s, only a single set of L exponential operations 

for the first value of k is needed. The other exponential terms can be derived from that of the 
previous value of k , with each derivation involving only -3 multiplication operations. Because a 
multiplication operation is much faster than an exponential operation, flux calculations can be 
greatly accelerated by choosing values of k's according to (4.12) for computing transmittances. 
The first absorption coefficient, k h the constant T|, and the k -distribution function are given in 

Table 4 for the water vapor. 


Table 4. Parameters for the transmittance given by (4.10) due to water vapor line absorption. 
k { / |1 is the first absorption coefficient, and T| is the constant given in (4.12). A g is the k - 
distribution function (4. 1 1 ). Band 3 is treated separately in Section 4.5. Units of k are g 1 cm 2 . 



Band 1 

Band 2 

Band 4 

Band 5 

Band 6 

Band 7 

Band 8 

Band 9 

Band 10 


2.96e+l 

4.17e-l 

1 5.25e-4 

5.25e-4 

9.37e-3 

4.72e-2 

1.32e+0 

5.25e-4 

1.06e-l 


6 

6 

6 

6 

8 

9 

6 

16 

8 

Ag/ 

.2747 

.1521 

.4654 

.5543 

.5955 

.1958 

.0740 

.1437 

.3153 

Ag 2 

.2717 

.3974 

.2991 

.2723 

.2693 

.3469 

.1636 

.2197 

.4604 

Ag. t 

.2752 

.1778 

.1343 

.1131 

.0953 

.3147 

.4174 

.3185 

.1326 

Ag., 

.1177 

.1826 

.0646 

.0443 

.0335 

.1013 

.1783 

.2351 

.0798 

Agj 

.0352 

.0374 

.0226 

.0160 

.0064 

.0365 

.1101 

.0647 

.0119 

Ag 6 

.0255 

.0527 

.0140 

.0000 

.0000 

.0048 

.0566 

.0183 

.0000 
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MINOR ABSORPTION BAND 


In addition to the major absorption due to water vapor, C0 2 and 0 3 , there are also many 
weak absorption bands of N 2 0, CH 4 , C0 2 , and CFCs located mostly in the water vapor window 
region between 800 and 1300 cm Individually, the effect of these bands on fluxes is small, but 
collectively the effect is to reduce the outgoing LW radiation at the top of a clear atmosphere by 
~4 W m 2 and to enhance the surface heating by ~2 W m 2 . These effects are significant in the 
study of anthropogenic climate change. 

The spectral bands which include the minor absorption are shown in Table 1. A total of 13 
sets of minor transmission functions are computed: four in Band 4, two in Band 5, four in Band 6, 
two in Band 7, and one in Band 10 (Sub-band 3a), 

The ^-distribution method is also applied to computing the transmission functions averaged 
over individual minor spectral bands (Table 1). The reference pressure and temperature, values of 
the parameters used for the minor bands are listed in Tables 5-7 (see Kratz et al., 1993, and Kratz 
et al., 1998, for details of the parameterization). Because the absorption is weak in those minor 
bands, it requires only a small number (1-4) of the /c-intervals for accurate transmission 
calculations. 
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Table 5. Parameters for computing the flux transmission functions due to N 2 0 and CH 4 
using the ^.-distribution method. (p r = 500 hPa and 0 r = 250 K). 




n 2 o 


ch 4 



Band 6 

Band 7 

Band 10 

Band 6 

Band 7 

A| / M 

6.32e-2 

5.36e-2 

2.52e-l 

5.81e-3 

6.29e-2 

n 

21 

8 

58 


12 

As, 

.9404 

.5620 

.971 

1.000 

.6107 

A*; 

.0596 

.1387 

.029 


.2802 

Air. 


.2406 



.1073 



.0587 



.0018 

nt 

.0 

.48 

.0 

.0 

.65 

a 

1 ,93e-3 

1.38C-3 

1.45e-3 

1.70e-2 

5.96e-4 

P 

4.38e-6 

7.48e-6 

3.67e-6 

1.58e-4 

-2.29e-6 


Table 6. Same as Table 5, except for the absorption due to C0 2 . 



Band 4 

Band 5 

Band 10* 

*l/|* 

1 .92e-7 

1 ,92e-7 

2.66e-5 

T\ 

5 

5 

8 

Agi 

.1216 

.0687 

.267 

2 

.2436 

.1480 

.220 

Agj 

.2498 

.1951 

.210 

Aj?4 

.2622 

.3344 

.241 


.0781 

.1720 

.020 

Ag* 

.0427 

.0818 

.042 

m 

.0 

.0 

.5 

a 

3.58e-2 

3.43e-2 

1.79e-2 

p 

4.05e-4 

3.74e-4 

1.02e-4 


* p r = 300 hPa, 
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Table 7. Same as Table 5, except for the minor CFC absorption bands. There is only one 
^-interval for each CFC band, and there is no pressure scaling applied, i.e. m=0. 



CFC- 11 

CFC-12 

CFC-22 

Band 4 

Band 5 

Band 4 

Band 6 

Band 4 

Band 6 

k 

1.90e+l 

1.02e+l 

1.58e+l 

3.70e+l 

6. 1 8e+0 

3.28e+l 

F 







a 

1.27e-3 

8.19e-4 

8.77e-4 

8.62e-4 

9.65e-4 

-3.00e-5 

P 

3.56e-6 

4.68e-6 

-5.88e-6 

-4.23e-6 

1.31 e-5 

5.25e-7 


4.2 Pre-computed Transmission Tables 

Although the ^-distribution method with linear pressure scaling is computationally very 
fast, it is not accurate in the middle atmosphere (0.01-10 hPa), where the pressure ranges by 3 
orders of magnitude and where the Doppler broadening of absorption lines is important. Cooling 
in the middle atmosphere is primarily due to C0 2 in the 15 -pm band (Band 3), in the 9.6 -pm 

band (Band 5) and secondarily due to water vapor near the centers of absorption bands (Bands 1 , 
2, and 8). 

It is shown in Chou and Kouvaris (1991) and Chou and Suarez (1994) that transmittances 
in these bands can be simply derived from pre-computed transmission tables. The basis for this 
method is that a nonhomogenous layer with pressure and temperature varying with height can be 
treated as if it were homogeneous with an effective pressure and temperature defined by 



(4.13) 


(4.14) 
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where u is the absorber amount, and the integration is over the depth of a layer. It has been well 
recognized (e.g. Wu, 1980; Chou and Kouvaris, 1991) that the contribution to the atmospheric 
cooling is primarily from near-by layers. The pressure and temperature variations among near-by 
layers are small, and the simple scaling approximations of (4.13) and (4.14) can be applied to 
accurate calculations of cooling rate. 


With the two-parameter scaling of (4.13) and (4.14), the flux transmission becomes a 
function of the absorber amount and the effective temperature and pressure. These dependencies 
can be accurately pre-computed from the following equation using a line-by-line method, 


T{u, Peff $eff ) 



(u,PeJf & e ff )Bv (0 o )dv 

f 


(4.15) 


where Av is the width of the entire spectral band. The band transmission varies rapidly with 

pressure, but rather smoothly with temperature. The size of the three-dimensional transmission 
tables can be reduced to three two-dimensional tables using a quadratic fit in temperature, 

T(u,p e ff,Q e ff) = a(u,p e ff) + b(u,p e ff)(Q e ff - 250) + c(w, p eff )(d eff -250) 2 (4.16) 

where a, b, and c are regression coefficients. This regression is valid for temperatures ranging 
from 170AT to 330 K\ for this range it introduces only ~ 1% error in the absorptance. Table 8 
shows the first values, the intervals, and the interval numbers of log w (u) and log l0 (p eff ) of the 
tables a , b and c for transmissions due to C0 2 absorption in the 15 pm region, 0 3 absorption 
in the 9.6 pm region, and water vapor absorption in the three strong absorption bands. It is 
noticed that the transmission tables can be applied to an atmosphere with a C0 2 concentration 100 
times of the present value, i.e. 35000 ppmv. 

The absorption due to 0 3 covers a spectral range wider than that of Band 5 shown in Table 
1. To include the weak absorption outside this spectral band, we have made adjustment to the 
transmission function by requiring that 

_ f v 2 r v t> 

(1-T) ^(7Vv= 4 -T v (T)dv (4.17) 

Jv i 
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Table 8. The first values, logjJw^j and log w (p eff ) h the intervals, Alog 1( Jw eff ) and 
Alog }0 (p eff ), and the numbers, nw and np , of log w (w eff ) and logJp eff ) of the tables a , b 

and c for transmittance calculations using (4.16). Units of w eff are g cm for water 
vapor and (cm-atm) STP for C0 2 and 0 3 , and the unit of p eff is hPa. 


Band 1 Band 2 Band 3 Band 5 Band 8 


Absorber 

h 2 o 

logiJW'ff), 

-8 


-2 

MogJw et f) 

0.3 

Alog,„(p rff ) 

0.2 

nw 

31 

np 

26 


H 2 0 

C0 2 

-8 

-4 

-2 

-2 

0.3 

0.3 

0.2 

0.2 

31 

30 

26 

26 


O 3 

H 2 0 

-6 

-8 

-2 

.9 

0.3 

0.3 

0.2 

0.2 

21 

31 

26 

26 


where v, - v 2 is the spectral range of Band 5, v a - v b is a wider spectral range that covers all the 
molecular lines in the absorption band, and x is the effective mean transmission of Band 5. Thus, 
x includes the effect of the 0 3 absorption outside the spectral range v a - v b . We could have 
expanded the range v, - v 2 , but the accuracy of the transmission calculations would be reduced 
when (4.27) is applied to include the water vapor effect, which is discussed in Section 4.4. The 
adjusted mean transmission is then pre-calculated from 

I (l-x v )By (T)dv 

; = (4.18) 

r v 2 

~R,(T)dv 

The method of table look-up is simple and accurate. As shown in Chou and Kouvaris 
(1991) and Chou and Suarez (1994), it can be applied to accurate calculating the fluxes and 
cooling rate in both the middle and lower atmosphere extending from 0.01 hPa to the earth's 
surface. However, this method could be significantly slower than the ^-distribution method on 
computers that cannot efficiently perform table look-ups. We have therefore provided the option 
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of using the ^-distribution method for all bands, except for the computation of O * transmission in 
Band 5. 


43 One-parameter scaling for water vapor continuum absorption 


As shown in Table 1, water vapor continuum absorption is included in Bands 3-7. The 
water vapor continuum absorption coefficient increases with increasing water vapor partial 
pressure but with decreasing temperature. It can be approximated by 





\Po) 




(4.19) 


where p e is the water vapor partial pressure in hPa , p () = 1013 hPa , and k vo is the absorption 
coefficient when p e = 1013 hPa and 0 = 296 AT. Values of k v0 are computed from the analytical 

representation given by Roberts et al. (1976), which is a fit to the laboratory data of D. E. Burch 
for Bands 3-6 and from Clough et al. (1989) for Band 7. 

With the scaling of (4.19) for the absorption coefficient, the flux transmission reduces to 


Tv(p,p')= 2 e 


t 


-K.o wC (P-P">k 




(4.20) 


where w‘ is the scaled water vapor amount for continuum absorption given by 


H 1 (p,p 


, ) = L 61 rp(£' 

S Jp\Poj 


i8oq 




] 


i 


6 (p") 296 


dp" 


(4.21) 


Here we have used the approximation 


dP 

Pe ~ 

F 0.622 

and q is the specific humidity. 

Taking the Planck-weighted average of (4.20) over a band we have 


J 1 

;K) 

1 £y(6 0 )^v 

J 

ft 

B, (Q a )dv 

Av 


(4.22) 


(4.23) 
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where we note that the transmission is a function of w c only. The absorption coefficient k vo varies 
slowly with wavenumber, except in Band 3 where it varies by a factor of =3.5. Therefore, we 
divide Band 3 into three sub-bands (see Section 4.5 for more detail of the division of Band 3 into 
3 subbands). The Planck-weighted flux transmissions for the three sub-bands and for Bands 4 - 7 
were computed from (4.23). These transmittances were then fit by 

T(w c ) =e -iI H ' t/M ( 4.24 ) 

The effective absorption coefficients, lc c / \L, for water vapor continuum absorption are given in 
Table 9. 


Table 9. The effective absorption coefficient lc c / fi for the water vapor continuum 
absorption for Bands 3-7. Units of k c / p. are g 1 cm 2 . 



Band 3a 

Band 3b 

Band 3c 

Band 4 

Band 5 

Band 6 

Band 7 

£ 

4 

109.6 

54.8 

27.4 

15.8 

9.40 

7.75 

8.78 


4.4 Overlapping of Gaseous Absorptions 

When there are more than one absorber involved in a spectral band, overlaps must be 
considered. The total transmission for, say, two absorbers, can be written as 

\T j(v)T 2 (v 

T n = l (4.25) 

r 

where the subscripts 1 and 2 denote the two absorbers. If the transmittance is expressed as the 
sum of the mean over the spectral band, T , and the deviation from the mean, T, then (4.25) 
reduces to 
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(4.26) 


Tn = 


f(ziH-r' 1 (v)X7' 2 +7’ i 2 (v))^ __ frjWF^v^v 

w <T» rrr J 


J 


= T\ T 2 + : 


dv 


J 


dv 


If the overall shapes of the absorption curves due to both absorbers are uncorrelated, the 
last term of (4.26) can be neglected, and the total transmittance becomes 

T\ 2 = T\T 2 (4.27) 

4.5 Special Treatment of the 15- pm Band 

The 15 -pm band poses additional difficulties for two reasons. First, the water vapor line 
absorption and the continuum absorption are highly correlated. As can be seen in the bottom two 
panels of Figure 2, the absorption coefficient increases with decreasing wavenumber by a factor 
of about 100 for water vapor line absorption and about 10 for continuum absorption. Thus, the 
water vapor line and continuum absorption are highly correlated, and the approximation (4.27) 
cannot be applied directly to the entire band. Second, the C0 2 absorption coefficients differ by 
several orders of magnitude between the band center and the wings (see the top panel of Figure 
2). Rather than trying to parameterize the correlation effect or the variations in C0 2 absorption, 
we simply divide the band into three sub-bands (see the bottom panel of Figure 2 and Table 1 ) 
and then combine the parameterized transmittances of the sub-bands into a single band 
transmittance. This transmittance is then used in the usual way to solve the transfer equations for 
the entire band. 
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c a 

2 




Wavenumber (cm 1 ) 


Figure 2. Schematic of the absorption coefficient for C0 2 and for water vapor line and 
continuum absorption in the 15 -fjm region (Band 3). The subdivision of the band into three sub- 
bands is made to accurately account for the large variations of the absorption coefficients. 


The H 2 0 line absorption in Band 3 is computed from the ^-distribution method. We divide 
the band into three sub-bands. Within each of the three sub-bands, the transmittances due to line 
and continuum water vapor absorption are sufficiently uncorrelated that we can overlap them 
using the multiplication approximation (4.27). Letting T, L and T, r be the line and continuum 

transmittances of sub-band /, we obtain the total water vapor transmittance in Band 3 by 
overlapping them and taking the Planck-weighted average 


jWV 



Ti L Ti C 1 


Wo). 


(4.28) 


where B t (Q 0 ) is the Planck flux integrated over Sub-band i, and B(Q 0 ) is the Planck flux integrated 
over the three sub-bands. 
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Substituting for the line transmittance from (4.10) and for the continuum transmittance 
from (4.24), we have 


T »v^ e -k^r^ e - K «r f( 

i - 1 n = 1 I 

i— 1 n =1 


Me) 

Wo) _ 


(4.29) 


We use a single set of k h and r\ for the three water vapor sub-bands. See (4.12) for relating kj to 

k n . Table 10 shows the values of these parameters. The pressure and temperature scaling 
parameters are shown in Table 3. Since we use the same scaling parameters for the three sub- 
bands, the scaled water vapor amounts, vv and w c , are independent of sub-band. Therefore, the 
exponentials for the water vapor line transmittance in the inner summation of (4.29) only need to 
be evaluated once. 


Table 10. The first value of the absorption coefficient / fi for the water vapor line 
absorption, the value of T| that k n = l\k n h and the Planck-weighted k -distribution 
function (Ag) given by (4.29) and (4,30) of the three sub-bands in the 15 fum region. 



WATER VAPOR 


CO, 

Subband a 

Subband b 

Subband c 

Wings 

Center 

k\i\y 

1.33e-2 

1 ,33e-2 

1.33e-2 

2.66e-5 

2.66e-3 


8 

8 

8 

8 

8 

(Ag)i 

.1782 

.0923 

.0000 

.1395 

.0766 

(Ag) 2 

.0593 

.1675 

.1083 

.1407 

.1372 

(Ag)3 

.0215 

.0923 

.1581 

.1549 

.1189 

(Ag) 4 

.0068 

.0187 

.0455 

.1357 

.0335 

(A g\ 

.0022 

.0178 

.0274 

.0182 

.0169 

( a «) 6 

.0000 

.0000 

.0041 

.0220 

.0059 


The continuum absorption coefficients, k c \ are given in Table 9. Note that the coefficients 
for Sub-bands 3a and 3b are integer multiples of the coefficient for Sub-band 3c; thus only one 


22 






exponentiation for the water vapor continuum transmittance needs to be performed in evaluating 
T w for the three sub-bands. 

The C0 2 transmittance may be computed either using the /.-distribution method or using 
table look-up. When the table look-up is used, a single transmittance is computed directly for the 
entire band. When the jt-distribution method is used, we separate Sub-band b (the center region) 
from Sub-bands 3a and 3c (the wings). The optical properties of Sub-bands 3a and 3c are very 
similar, and so we have combined them by using the same scaling parameters (p n 0„oc, P, and m 

in Table 3 and k, in Table 10). The band-averaged C0 2 transmittance for Band 3 is thus computed 
from 


T c ° 2 <«°> 

n=l n = 1 

where the subscripts c and w denote, respectively, the band center (Sub-band 3b) and band wings 
(Sub-band 3a and 3c), and u c and u w are the scaled C0 2 amount in the band center and band 
wings. Note that 

^( A £c),,+( A sJ n = 1 (4 - 31) 

/!= 1 


Values of (Ag) are given in Table 10. 

Finally, the total flux transmittance in the 15-jim band is computed from 

ji rjiWV jrCOj 


(4.32) 


5. SPECIAL TREATMENT OF THE \l-fim N 2 0 BAND 

The spectral range of the absorption due to N 2 0 in the 11 -pm region is relatively narrow 
(560-620 cm 1 ) and coincides approximately with the Sub-band 3a. As shown in Figure 2, the 
absorption in this subband differs significantly from Sub-bands 3b and 3c. Therefore, the 
absorption due to N 2 0 in this band is treated differently from the others and is designated as Band 
10 shown in Table 1. Only the changes in fluxes due to the N 2 0 absorption are computed. From 
(3.9) and (3.10), the change in fluxes can be computed from 
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(5.1) 


-f 


AF3a (/»)= 5 3a (0’) 


MM£±}\ dp . 

op' 


A/ 7 3fl (/>) ®e/03fl(0*)A73 fl (/>,p J ) + 


J, 


dp 


(5.2) 


In the above equations, the change in the mean transmittance due to N 2 0 is defined as 

A73 a =Tla- T-ia = T L 3aT C 3aT C ° 2 3 a [T N i° 3 a - 1] (5.3) 

where the subscript 3a denotes Sub-band 3a, the superscript * denotes inclusion of the N 2 0 
absorption, and the superscripts L, C, C0 2 and N 2 0 denote water vapor line, water vapor 
continuum absorption, C0 2 and N 2 0 respectively. The upward and downward fluxes are then 
adjusted according to 

F* (p)= F^{p)+ AFiJ’(p) (5.4) 

F* (p) = f\p)+AF^(p) (5.5) 

The transmission functions are computed using the ^-distribution method (4.10) for T L , 

T C0 2 and ant j t h e one-parameter scaling method (4.24) for T C . The parameters used for 

computing the transmission functions of Band 10 are given in Tables 3-6. 


6. CLOUDS 

6.1 Cloud Single-Scattering Properties 

Scattering by clouds is parameterized separately for water droplets and ice particles. The 
high-spectral resolution extinction coefficient, the single-scattering albedo, and the asymmetry 
factor were computed from the Mie theory for water droplets and the method of Fu and Sun 
(1998) for hexagonal ice particles. Dr. Si-Chee Tsay of NASA/Goddard Space Right Center 
provided us data of these parameters for water drops, and Prof. Qiang Fu of the University of 
Washington provided us the data for ice particles. Mean values of the extinction coefficient (3, 

single-scattering albedo to, and asymmetry factor g are in turn derived for each spectral band 
according to 
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(6.1) 


P( * ) = ^ (3v ( r e ) B, (0 O ) / ^ B v (0 O ) 

Av Av 

0)( r e ) = ^co v (r e ) Pv, ( r e )Bj (0 O ) Pv (*)* (0 O ) (6-2) 

Av Av 

ft r ,. ) = ^ ft, (r e )co v (r f ) pv (r*)#, (0 O ) / ^co v (/- e )Pv (r e )^, (0 O ) (6.3) 

Av Av 

where r, is the effective size of cloud particles which is a function of the ratio of the particle 
volume to the surface area , B is the Planck function, 0„ is a medium value of the atmospheric 

temperature set to he 250 K, and Av is the spectral interval of a band. We fit the mean single- 
scattering parameters by the following functions of r e , 

C ft 2 

P< r , ) - cu | + — for ice particles (6.4a) 

^ < / \CR t 


= 2j c v jre 

;= l 

4 

(Mre)=^' j Cujr e 

i= 1 


i- 1 


i-l 


4 

8(r e )=^c g jr e l ~ X 
i = 1 


for liquid drops 


for both ice and liquid water particles 


for both ice and liquid water particles 


(6.4b) 


(6.5) 


( 6 . 6 ) 


Tables 11-13 list values of the regression coefficients, c ^ c g , for both water droplets and ice 
particles. Finally, the cloud optical thickness is derived from 

x= p CWP (6.7) 


where CWP is the cloud water path (column mass per unit area), either ice or liquid water. In the 
current version of the code, the effective radii r e are input parameters and are not computed 
within the code. 
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Table 1 la. Coefficients for computing the extinction coefficient for cloud 
ice particles from (6.4a). Units of the particle size in (6.4a) are pm. 


Band 

C P ./ 

C P2 


i 

-0.44171 

0.62951 

0.06465 

2 

-0.13727 

0.61291 

0.28962 

3 

-0.01878 

1.67680 

0.79080 

4 

-0.01896 

1.06510 

0.69493 

5 

-0.04788 

0.88178 

0.54492 

6 

-0.02265 

1.57390 

0.76161 

7 

-0.01038 

2.15640 

0.89045 

8 

-0.00450 

2.51370 

0.95989 

9 

-0.00044 

3.15050 

1.03750 

10 

-0.02956 

1 .44680 

0.71283 


Table 1 lb. Coefficients for computing the extinction coefficient for cloud 
liquid water droplets from (6.4b). Units of the particle size in (6.4b) are pm. 


Band 

C V 

C P2 


c B-/ 

i 

0.08641 

0.01769 

-1.5572E-3 

3.4896E-5 

2 

0.22027 

0.00997 

-1.8719E-3 

5.31 12E-5 

3 

0.38074 

-0.03027 

1.0154E-3 

-1.1 849E-5 

4 

0.15587 

0.00371 

-7.7705E-4 

2.0547E-5 

5 

0.05518 

0.04544 

-4.2067E-3 

1.0184E-4 

6 

0.12724 

0.04751 

-5.2037E-3 

1.371 IE-4 

7 

0.30390 

0.01656 

-3.527 IE-3 

1.0828E-4 

8 

0.63617 

-0.06287 

2.2350E-3 

-2.3177E-5 

9 

1.15470 

-0.19282 

1.2084E-2 

-2.5612E-4 

10 

0.34021 

-0.02805 

1.0654E-3 

-1.5443E-5 
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Table 12a. Coefficients for computing the single-scattering albedo for 
cloud ice particles from (6.5). Units of the particle size in (6.5) are pm. 


Band 


C <4 2 


C 0 

1 

0.17201 

1.2229E-2 

-1.4837E-4 

5.8020E-7 

2 

0.81470 

-2.7293E-3 

9.7816E-8 

5.7650E-8 

3 

0.54859 

-4.8273E-4 

5.4353E-6 

-1.5679E-8 

4 

0.39218 

4.1717E-3 

-4.8869E-5 

1.9144E-7 

5 

0.71773 

-3.3640E-3 

1.9713E-5 

-3.3189E-8 

6 

0.77345 

-5.5228E-3 

4.8379E-5 

-1.5151E-7 

7 

0.74975 

-5.6604E-3 

5.6475E-5 

-1.9664E-7 

8 

0.6901 1 

-4.5348E-3 

4.9322E-5 

-1.8255E-7 

9 

0.83963 

- 6.7253E-3 

6. 1 900E-5 

-2.0862E-7 

10 

0.64860 

-2.8692E-3 

2.7656E-5 

-8.9680E-8 


Table 12b. Coefficients for computing the single-scattering albedo for 
cloud liquid water droplets from (6.5). Units of the particle size in (6.5) 
are pm. 

Band 


C 0)2 



1 

0.17201 

1.2229E-2 

-1.4837E-4 

5.8020E-7 

2 

-0.01338 

9.3134E-2 

-6.049 IE-3 

1.3059E-4 

3 

0.03710 

7.321 IE-2 

-4.421 IE-3 

9.2448E-5 

4 

-0.00376 

9.3344E-2 

-5.6561E-3 

1.1387E-4 

5 

0.40212 

7.8083E-2 

-5.9583E-3 

1.2883E-4 

6 

0.57928 

5.9094E-2 

-5.4425E-3 

1.2725E-4 

7 

0.68974 

4.2334E-2 

-4.9469E-3 

1.2863E-4 

8 

0.80122 

9.4578E-3 

-2.8508E-3 

9.0078E-5 

9 

1.02340 

-2.6204E-2 

4.2552E-4 

3.2160E-6 

10 

0.05092 

7.5409E-2 

-4.7305E-3 

1.0121E-4 
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Table 13a. Coefficients for computing the asymmetry factor for cloud ice 
particles from (6.6). Units of the particle size in (6.6) are pm. 


Band 

c i.i 



C S ' 4 

1 

0.57867 

1.0135E-2 

-1.1 142E-4 

4.1537E-7 

2 

0.72259 

3.1 149E-3 

-1.9927E-5 

5.6024E-8 

3 

0.76109 

4.5449E-3 

-4.6199E-5 

1 .6446E-7 

4 

0.86934 

2.7474E-3 

-3.1301E-5 

1.1959E-7 

5 

0.89103 

1.8513E-3 

-1.6551E-5 

5.5193E-8 

6 

0.86325 

2.1408E-3 

-1.6846E-5 

4.9473E-8 

7 

0.85064 

2.5028E-3 

-2.0812E-5 

6.3427E-8 

8 

0.86945 

2.461 5E-3 

-2.3882E-5 

8.243 IE-8 

9 

0.80122 

3.1906E-3 

-2.4856E-5 

7.241 IE-8 

10 

0.73290 

4.8034E-3 

-4.4425E-5 

1.4839E-7 


Table 13b. Coefficients for computing the asymmetry factor for cloud liquid 
water droplets from (6.6). Units of the particle size in (6.6) are pm. 


Band 

C g.l 

c tJ 

C *J 

C g,4 

1 

-0.51930 

0.20290 

-1.1747E-2 

2.3868E-4 

2 

-0.22151 

0.19708 

-1.2462E-2 

2.6646E-4 

3 

0.14157 

0.14705 

-9.5802E-3 

2.0819E-4 

4 

0.41590 

0.10482 

-6.91 18E-3 

1.5115E-4 

5 

0.55338 

7.7016E-2 

-5.2218E-3 

1.1587E-4 

6 

0.61384 

6.4402E-2 

-4.624 IE-3 

1.0746E-4 

7 

0.67891 

4.8698E-2 

-3.702 IE-3 

9.1966E-5 

8 

0.78169 

2.0803E-2 

-1.4749E-3 

3.9362E-5 

9 

0.93218 

-3.3425E-2 

2.9632E-3 

-6.9362E-5 

10 

0.01649 

0.16561 

-1.0723E-2 

2.3220E-4 
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6.2 Layer optical properties 


When the absorption and scattering due to clouds interact with that of gases and aerosols in 
a layer, the effective optical thickness, single-scattering albedo, and asymmetry factor in a 
spectral band are computed from 


T= > T/ 


(0 = 




/ i 




( 6 . 8 ) 


(6.9) 


( 6 . 10 ) 


where the summation is over all gases, aerosols, and liquid water and ice particles, i. Thus, in the 
LW radiation code, ice and liquid water cloud particles are allowed to mix in a layer. 


6.3 Paramaterization for Cloud Scattering 

Explicit calculations of scattering by clouds requires a large amount of computer time. 
Instead, we scale the cloud optical thickness x to approximate the effect of scattering according 

to. 


x c = ( l— co fyz (6.11) 

where / is the fraction of radiation scattered downward (upward) for radiation incident from 
above (below). This scaling takes into account the effect of back-scattering and absorption by 
cloud particles. 

Assuming that the scattering phase function can be approximated by the Henyen- 
Greenstein function, the forward-scattering function /is computed as a function of the asymmetry 
factor and fit by a polynomial function 

(=1 

where a, = 0.5, a 2 = 0.3738, a } - 0.0076, and a 4 - 0.1 186. 


29 



Thus, the cloud optical thickness is reduced from x to x c , and scattering is not explicitly 
computed. By scaling the cloud optical thickness (6.11), the effect of scattering is parameterized 
but not explicitly computed. The computing time is the same as the case without taking into 
account the cloud scattering. See Chou et al. (1999) for details of the parameterization for the 
cloud LW scattering. 

6.4 Cloud Overlapping 

Equation (2.3) applies to the case where the atmosphere is horizontally homogeneous. 
When clouds with fractional areal cover occur at different heights, assumptions have to be made 
on the ways these clouds are overlapped for radiative flux calculations. There are a number of 
schemes used to vertically associate clouds at different height (e. g. Manabe and Strickler, 1964; 
Geleyn and Hollingsworth, 1979; Liang and Wang, 1997; Raisanen 1998; Li., 2000). These 
schemes all assume that clouds either are randomly overlapped or have a combination of random 
and maximum overlapping. In our radiation parameterization, we also adopt a maximum-random 
overlapping cloud scheme. 

If a cloud layer with fractional cover A and optical thickness x/ is introduced between p 
and p\ the mean radiance transmittance in the direction p becomes 

[l _ j\]e~ Zv + Ae~^ v (p - p )+z ' /C(p ’ p ' ] /fi 

and the mean flux transmittance is 


1 

T, * (p, p') = 2 jj [1 - A]e~ * (p - p ' > * + Ae~^ v (p ' p ' )+z A‘» 


(6.13) 


= 1 - 




T v (p,p') 


where x v is the optical thickness of a clear atmosphere, JH is the effective mean value of p which 
converts the radiance transmittance to flux transmittance, and T v is the transmittance of a clear 
atmosphere defined by (2.3). The diffusivity factor 1 /(I depends on x v but is commonly chosen 

to be 1.66. Letting M v -A(\-e~ Zy> ^ ), Equation (6.13) becomes 


Tv (pyp) — (l” N v )7^ (p>p') 


(6.14a) 
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If we further assume that x v c is constant within a spectral band and integrate Equation (6.14a) over 
the band, we have 

T*(p,p') = {\-N)T(p,p') (6.14b) 

where T is the band-averaged flux transmittance defined by (3.7). Because 1 - e is the flux 
emissivity of the cloud layer, N can be regarded as the effective cloud cover with an emissivity of 
1. Alternatively, (l-N) can be regarded as the effective cloud transmissivity with a fractional 
cloud cover of 1 . 

When there is more than one cloud layer with fractional cover between p and p\ the 
situation is considerably more complicated since we need to describe how the clouds are 
overlapped. In general we can write 

T*(p y p') = C(p,p')T(p,p') ( 6 - 15 ) 

where C(p,p') depends on the values of A and r c for the various cloud layers. 

For the special case of any combination of overcast and randomly overlapped fractional 
cloud layers, we have 

C(/>,/>')=]^[(l-A^) ( 6 - 16 ) 

j 

where the subscript j denotes the cloud layers between p and p’. Since n ( 1 — Nj) is the fraction 

of the horizontal area that would be cloud-free if all cloud layers (including overcast layers) were 
assigned their equivalent black-cloud fraction and then randomly overlapped. Harsh vardhan et al. 
(1987) referred to C(p,p') as the probability of a clear line-of-sight between p and p'. Obtaining 
C(p,p') for non-random arrangements of gray clouds, however, is not as straightforward. In 
general, C(p,p') does not depends on N alone, but on more complicated combinations of A and 

x c . 

Consider another arrangement discussed by Harshvardhan et al. (1987), that of maximally 
overlapped clouds. If the largest cloud between p and p' is black, then C(p,p') is simply equal to 
one minus its fraction. If all clouds are allowed to be gray, however, C(p,p') will depend on all 
cloud fractions and optical thicknesses. For the maximally overlapped clouds between p and p\ 
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the clear line-of-sight C(p,p *) is independent of the order of cloud layers and may be computed by 
putting the clouds in order of increasing cloud cover: 


For the special case with three maximally overlapped clouds between p and p\ the clear 
line-of-sight C(p,p ’) may be computed by putting the clouds in order of increasing cloud cover: 


1 - C = Ai 


1 —e 


+ (^2 ~ ) 


1 -e 


( -(X C 2+T C 3) 


= A 3 ( l-e~ x 3 )+ 


A 2 ( \-e x 2 ) + A\(\-e' 


+ (A 3 -A 2 ) I 

Tf| ) e ~ xC 2 Y z ‘ 3 


1 —e 


,-T 3 


(6.17) 


= (V-j + 


N 2 +N x e~ T 2 If 


-t 3 


where In (6.17), we have dropped \x from x c/ ^ for simplicity. As referenced to the upper 

cloud group in Figure 3, the first and second brackets represent the absorption of the regions 
where there are three and two cloud layers overlap, while the third bracket represents only the 
absorption of the largest cloud cover. 

It follows from (6.17) that for the case with J maximally overlapped clouds, C(p,p r ) may be 
computed by evaluating the following recursion: 


D {0) = 0 



(6.18) 

D (J) = Nj + D (j ~ l) e~ zCj,ii , 

j = 1, ■ 


(6.19) 


and 

C(p,p') = 1-D {J) (6.20) 

Whenever a black cloud (x c j— >oo) occurs in (6.19), all smaller clouds are eliminated from the 

recursion. This treatment of cloud overlapping was first proposed by Chou and Suarez (1994). A 
similar cloud overlapping scheme is also proposed by Raisanen (1998). 

To provide an example for the overlapping scheme, we show in Figure 3 two groups of 
maximally overlapped clouds between p and p\ The upper group and the lower group are labeled 
U and L, respectively. These two groups are assumed randomly overlapped. The clear line-of- 
sight between p and p’ is 

C(p,p')=C u C 1 (6.21) 
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where C' 7 and C* are the dear line-of-sight of the upper and lower cloud groups computed from 
(6.20). Fluxes can then be computed from (3.9) and (3.10) with T(p.p') replaced by C(p,p')T(p,p'). 

In the current version of our radiation code, clouds are grouped into three height ranges: 
high, middle and low separated approximately by 400 hPa and 700 hPa pressure levels. Clouds 
within each height group are assumed maximally overlapped, and clouds among the three height 
groups are assumed randomly overlapped. This cloud overlapping scheme can be easily extended 
to include a more flexible situation that assumes maximum overlapping of adjacent cloud layers 
and random overlapping of cloud layers separated by clear layers, as originally proposed by 
Gelygn and Hollingswoth (1979). 



c 


u 


C L 


Figure 3. Schematic showing two groups of randomly overlapped clouds. Cloud layers within each group 
are assumed maximally overlapped. C v and C L are the clear line-of-sight of the upper group and lower 
group, respectively. The clear line-of-sight between the levels p and /?' is C V C L as shown in (6.21 ). 


1 . AEROSOLS 

The aerosol optical thickness, single-scattering albedo, and asymmetry factor for each of 
the 10 bands and each of the atmospheric layers are specified input parameters. Different types of 
aerosols are allowed to co-exist in a layer. The effective optical thickness, single- scattering 
abledo, and asymmetry factor of a layer are computed similarly to that of clouds, 
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(7.1) 


i 

i i 

g a = y.g fl ia> a »t a //y i J <o a /t fl » (7-3) 

f i 

where i is the index for aerosol type. The same parameterization for cloud scattering, (6.1 1) and 
(6.12), is applied to an aerosol-laden layer. We further assume that the aerosol layer is plane- 
parallel and the flux transmission function of the layer is computed from 

T a = T e x “ ^ (7.4) 

where x* a is the scaled aerosol optical thickness according to (6.11) and (6.12), and T and T are 
flux transmittances with and without aerosols. 


8. VERTICAL DISCRETIZATION 


8.1 Downward and Upward Fluxes 


To approximate the vertical integrals in (3.9) and (3.10), the atmosphere is divided into L 
layers numbered as shown in Figure 4. The downward and upward fluxes at level / for the ith 
band are computed from 


/-I 

F U l = Y B i[Ti(l,r+\)-Ti(l,n] 

^ id -H“ 

/=1 2 

L 

ij’4— 

r=i 2 

+ T i (l,L+\)F iL+ \ 


1 = 2, L+l ( 8 . 1 ) 


/ = 1 , L (8.2) 
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1 


2 


> 

^ F/ B/- 1/2 


Bl + l / 2 

t 


/-I 

/ 

/+ 1 


T i. r 



L 



Figure 4. The vertical grid and placement of various quantities for an atmosphere consisting of L layers. 
Quantities defined at the layers, such as the Planck flux, B are denoted by half-integer subscripts, and 
quantities defined at the levels separating layers, such as the downward and upward fluxes F, by integer 
subscripts. The transmittance shown is for a multi-layer region bounded by levels / and /'. Note that the 
surface is treated as a fictitious layer at L+3/2. 
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(8.3) 


Fi/ = 0 

^ix + l t =e i B i (B s ) + {l-E i )F u+ i i (8.4) 

where B u+l/2 is the Planck flux at the temperature 0 ;+//2 integrated over the band /, 7’//,/') is the 

Planck- weighted flux transmittance between the levels / and T,(l,L+l) is the transmittance 
between the level / and the Earth's surface, and 0, is the Earth's surface temperature. 


By rearranging terms in Equations (8.1 ) and (8.2) and defining 


AF+(lJ) = Ti(l,r) 




B , I"*-, 1 

/,/ — i ,l+- , 

2 2 ) 


r>i, 1=1,2 l 


( 


Af} T (/,/■) = Tj(l,F) 




B r '~ B r 1 

V 2 2 y 




B i =0 
/ — 

2 


3 

i X+"" 
2 


the downward and upward fluxes can be written as 

/-l 




2 /=i 


/=2, 3, ... , L+l 


L+l 


fJ = B , +7;'(/,L + l)(l-e 1 )^x+i' L - y A/- T (/,/') 1=1,2,... ,L 


/*=/+! 


(8.5) 

( 8 . 6 ) 

(8.7) 

( 8 . 8 ) 

(8.9) 

( 8 . 10 ) 


As opposed to (8.1) and (8.2), only one transmittance, T/l,l'), appears in each term under 
the summation in (8.9) and (8.10). The sums can be build-up term by term. When the 
transmittance of a layer bounded by the upper level l and the lower level l is computed, the 
downward flux at i and upward flux at / are immediately undated. In this way, there is no need to 
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store the transmittance matrix. The storage of the entire routine scales then like L , rather than Is 
— a very significant advantage for models with high vertical resolution. 

8.2 Treatment of Adjacent Layers 

Accurate calculations of the contribution of a layer to the flux at the boundaries of the layer 
arc very' important for two of reasons, (a) Among the atmospheric layers, the adjacent layers have 
the largest contribution to the fluxes, (b) The transmittance decreases rapidly within the adjacent 
layer. Significant errors might be induced if the emitting temperature of the adjacent layer is 
assumed to be the mean layer temperature, as in (8.1) and (8.2). To reduce the error arising from 
this assumption, we compute the emitting temperature (or Planck function) of the layer as a 
function of the layer transmittance and the temperatures at the upper and lower boundaries of the 
layer. 

The contribution to the upward flux at the level p from the adjacent layer below can be 
written as 

/• /?+ £p 

AF T (p)=-l B(p')[dT(p,p') /dp']dp' 

J p (8.11) 

= B U ( \-T 0 ) 

where A p is the thickness of the adjacent layer, B u is the effective Planck flux of the adjacent 
layer, and T„ = T(p + Ap, p) is the flux transmittance between p and p + Ap. 

The Planck flux B(p' ) varies slowly within a layer, but the transmittance T(p, pj 

decreases exponentially as p 1 departs from p. For T a — > 0, the adjacent layer is rather opaque, and 

the radiation primarily comes from regions close to the level p. The effective Planck flux B“ for 
this case should be close to B(p). On the other hand, if the adjacent layer is rather transparent (T 0 
— > 1), the transmitttance varies slowly with p' and the flux at p is contributed more or less evenly 

from the entire layer. By linearizing the Plank flux within the layer, it follows that B u — » B(p + 

Ap/2 ) for T„ -> 1. Calculations of the flux contribution from the adjacent layer should take into 

account the dependence of B u on T„. 

Assuming that the Planck function varies linearly with pressure and the Planck-weighted 
transmittance decreases exponentially with pressure away form p, we have B(p') = B(p) + 
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[B(p+Ap) - B(p)] (p'-p) / A p and T(p,p )-e cl p p \ where c is a constant. It follows from 
(8.1 1 ) that the effective Planck flux of the adjacent layer is 




p) — B{ p+ Ap)T n 

l-T„ 


B(p)-B(p+ Ap) 


lnT„ 


( 8 . 12 ) 


Equation (8.12) meets the conditions that B u -^B(p) for T 0 — >0, and =[ 6 (p + 4 ?)+ B(p)]/ 2 


for weak absorption when * cAp+ ^(cAp) 


Similarly, for the downward flux at the lower boundary of the layer, the effective Planck 
function of the layer is 


B(p+ Ap)-B{p)T n 


B(p+ Ap)-B(p) 

\-T 

L 1 l o J 

- 1 - 

lnT 0 


From (8.12) and (8. 1 3), we have 


(8.13) 


B^ + B u = B(p)+ B(p+ Ap) 

By replacing B u -. n f and B u+os with B J and B\ Equations (8.9) and (8.10) reduce to 

^ 1-2 


F i / = B d i,r-l+T i (rj-lj 


«... 1 

, lj— lj 

V 2 2 , 


H- 


]^A Ft {IX 


/ = 1 


f 


F[j - B u ,-7;.(/,/ + l)| 
ij f- 
2 




i-*, 3 

ij + — ij+- 
V 2 2 j 

L + 1 


+7j(/,L+l)(l-e J )fvx+l^ - ^ A Fi'dJ) 


1 = 1+2 


where 


g ‘*jr 


'Bu-BU'-iTidW-f 


Bi.r 

1-7; (/',/’-!) 

+ 

InT^FJ- 1) 


*“., =(S (i ,+B li ,.)-B'' , 


and /'=/+;. 


(8.14) 


(8.15) 


(8.16) 


(8.17) 

(8.18) 
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When a model layer spans a region where the temperature lapse rate changes signs, we 
cannot assume that B linearly varies with p, and Equations (8.17) and (8.18) are no longer valid. 
In such a case, we simply compute 3 d and B u from 


B u x =B d 
ij +— 



0.5 B |+ 0.25 (Bij+Bif) 
ij ■+— 


( 8 . 19 ) 


8.3 The k-Distribution Transmission Functions 

For gaseous absorption using the ^-distribution method, the transmittance of each of the k- 
intervals can be first computed for the L layers from 

AJ | = e ~ k n Aw i+ml^ 1=1,2 L; n=l,2, ...,N (8.20) 

/?,/+— 

2 

where k n is the absorption coefficient of the n' h ^-interval, and A vr ltK is the scaled absorber 

amount of the layer between the levels / and /+1. Both k and Aw vary with the spectral interval i. 

This index is dropped for simplicity. Because the set of k's are chosen according to (4.12), we 
need only to compute the exponential for the first value of k, and the others can be derived from 

Ar fl =(A7-„_i) n n = 2, 3, ..., N (8.21) 

Once AT, is computed, the computation of each of the A T„ with n> 1 requires only a few 
multiplication operations. 

The transmittance between an upper level / and an lower level /', i.e. I ’> /, is given by 


f-1 

T n (l,r) = J~\jAT i 1=1 T; T =1+1 L+l 

J. 1 n,j+- 

j=‘ 2 


and the transmittance between the level / and the surface is given by 


( 8 . 22 ) 


L 

7’„(/,L + l) = "rT AT , (8-23) 

X J. fl-f+T 

j=l 2 


Following (4. 10), the band-averaged transmission terms in (8.5) and (8.6) are computed from 
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( 8 . 24 ) 


N 

nu') = ^T n an Ag n 

n= 1 

8.4 Table Look-Up for Transmission Functions 

When using the pre -computed transmittance tables, the transmittance between the upper 
level / and the lower level V is derived from 

Tff = T(U(j.p{j\Qif) ( 8 . 25 ) 


where u is the (nnn-scaled) absorber amount, and p and 0 are the effective pressure and 
temperature. They are defined by 

/-I 




x 


Au i 
j+ - 


( 8 . 26 ) 


PIT 


r - i 


i A u 


H 


J+ 2 


Zj A u . 1 
4^ y-t— 

j=l 2 


6 /./' 



,A« , 

+— l+T 
2 J 2 


/-I 



( 8 . 27 ) 


( 8 . 28 ) 


where j+1/2 denotes the layer between the levels j and j+J. In (8.25), the transmittance T varies 
with the spectral band i. This index is dropped for simplicity. 
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9. SUB-GRID VARIATIONS OF LAND SURFACE 


If the size of a model grid-box is greater than tens of kilometers, the variations of surface 
temperature and emissivity within the grid box could have a significant effect on model 
simulations. In such a case, the emissivity and reflectivity of the model surface should include the 
effect of the sub-grid variations. Let us consider a vegetation-covered land surface characterized 
by a ground-skin temperature 0 g , vegetation canopy temperature 0,, ground emissivity e g , canopy 

emissivity £,., and canopy reflectivity r v . The radiation emitted by the ground and the vegetation 

that leaves the top of the canopy is given by 


T v \e g B(8 g ) + r^ v B(e v ) 

E = — ^ 2 J 

o-vp 


+e v B(0v) 


(9.1) 


where r g = (l-e g ) is the ground reflectivity, and T v = (/-£,.-r v ) is the canopy transmissivity. The 

first term within the bracket is the radiation emitted by the ground, the second term within the 
bracket is the radiation emitted downward by the vegetation and subsequently reflected by the 
ground, and the last term is the upward emission by the vegetation. The denominator in (9.1) 
takes care of the multiple reflections between the vegetation and the ground. 

The downward flux from the atmosphere is reflected by the vegetation and the ground. This 
term is given by 


R=rF i {p s ) 


(9.2) 


where r is the effective reflectivity of the canopy and the ground given by 


r- - r v +- 


(1- r v r g ) 


(9.3) 


For the case without vegetation cover, we have T v = 1 . E v = r v = 0, and Equations (9.1) and 
(9.3) reduce to 

£ = e g fl(0 g ) (9.4) 

r= r g (9-5) 

It should be noted that Equations (3.10), (8.2) and (8.4) are derived based on (9.4) and (9.5). 
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If there are N sub-grid boxes in a model box with fractional cover f n ,n= 1, 2, ... N, 
Equations (9. 1 ) and (9.3) become 


N 


■= 5 /" 






n= 1 


(1 r v,n r gjt) 




(9.6) 


N 

r= X /n 

«=1 

where 

N 

^ fn = 1 (9.8) 

n=l 

and the subscript n denotes the sub-grid box. 

For the case without vegetation cover, Equations (9.6) and (9.7) reduce to 
N 

E = fn£ g,nB(®g,n ) (9.9) 

n= 1 
N 

r= ^^^fn r g,n (9.10) 

n =1 

With E and r given above and R given by (9.2), Equations (8.2) and (8.4) reduce to the following 
general forms 

L 

f/ T = y B AniJ) -T{IS+\)]+T{1,L + \)F l ^ (9.11) 

/'+- 

/’=/ 2 

F L+l T = E+R=E+rF L+] l (9.12) 

Finally, the partial derivative of the upward flux with respect to the surface temperature given by 
(3.12) reduces to the following general forms 
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where B’(6) is the derivative of B with respect to temperature at 0. For simplicity, we have 
omitted the index for the spectral band i in all of the equations. 
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10. COMPARISONS WITH LINE-BY-LINE CALCULATIONS 


Fluxes and cooling rate computed using the transmittance parameterizations are compared 
with high spectral-resolution line-by-line calculations. Our line-by-line calculations of the 
absorption coefficient use the Air Force Geophysical Laboratory 1996 edition of the molecular 
absorption parameters (Rothman et al., 1998). The molecular line shape is assumed to follow the 
Voigt function. The absorption coefficient at wavenumbers > 10 cm 1 from the line center is 
taken to be zero, which is equivalent to a line-cutoff of 10 cm 1 . The absorption coefficient is 
computed at spectral intervals of 0.001 cm 7 for Band 3 and 0.002 cm 1 for the other bands. 

The mid-latitude summer atmosphere and the sub-arctic winter atmosphere taken from 
McClatchey et al. (1972) are used in the flux calculations. No clouds and aerosols are included. 
The C0 2 concentration is fixed at 350 ppmv , and the specific humidity above the tropopause is set 
to 4x1 0~ 6 g/g. The atmosphere is divided into 75 layers with A p ~ 25 hPa at pressures > 100 hPa 

and A log l(j p = 0.15 at pressures < 100 hPa. 

Cooling rates computed using the line-by-line method (solid lines) and the parameterization 
with the HIGH option (dashed lines) are shown in Figures 5 and 6. The cooling rate is presented 
as a function of log l0 p in the upper panels and p in the lower panels. It is noted that for the HIGH 
option, the transmission functions of C0 2 in Band 3 and of water vapor in Bands 1, 2, and 8 are 
computed using table look-up (see Table 1). It can be seen in the figures that, compared to the 
line-by-line method, the cooling rate is computed accurately. The maximum error in the 
stratosphere is < 0.25 °C/day and is < 0.15 °C/day in the troposphere. These errors are small when 
compared with the maximum cooling of ~12 °C/day in the stratosphere and ~TC/day in the 
troposphere. 

The downward flux at the surface and the upward flux at the top of the atmosphere are 
shown in Tables 14 and 15 for the mid-latitude summer atmosphere and the sub-arctic winter 
atmosphere, respectively. Errors are generally smaller than 1 Wm 2 in individual bands. 
Integrated over the LW spectrum, the parameterizations can compute the fluxes to within 1% of 
the line-by-line calculations for the mid-latitude summer atmosphere and within 0.5% for the sub- 
arctic winter atmosphere. 
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Cooling Rate (C/day) 


Figure 5. The cooling rate computed for a clear mid-latitude summer atmosphere using the line-by-line 
method (solid curves) and the "HIGH" option of the parameterization (dashed curves). Cooling is due to 
water vapor molecular line and continuum absorption, as well as C0 2 and (^absorption. 
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Table 14. Downward fluxes at the surface, F i s)n and upward fluxes at the top, , op , for 
a clear mid-latitude summer atmosphere computed using a line-by-line method and the 
parameterization with the HIGH option. Units of the fluxes are Wm \ 


Spectral Band {cm 1 ) 

F\< 

F r 

1 top 

Absorber 


0- 340 

line-by-line 

51.02 

34.56 

H 2 0 (line) 


Parameterization 

51.04 

34.40 

H 2 0 (line) 


340- 540 

line-by-line 

80.76 

61.01 

H 2 0 (line) 


Parameterization 

81.23 

60.54 

H 2 0 {line) 


540- 800 

line-by-line 

105.85 

67.90 

H 2 0 {line, continuum ), 

co 2 

Parameterization 

107.29 

68.29 

H 2 0 {line, continuum ), 

co 2 

800- 980 

line-by-line 

26.89 

58.62 

H 2 0 {line, continuum ), 

o 3 

Parameterization 

27.01 

58.72 

H 2 0 (line, continuum) 


980-1100 

line-by-line 

12.36 

22.36 

H 2 0 {line, continuum), 

o 3 

Parameterization* 

12.54 

21.41 

H 2 Q {line, continuum), 

Oj 

1100-1215 

line-by-line 

9.73 

20.96 

H 2 0 {line, continuum), 

o 3 

Parameterization 

9.63 

21.42 

H 2 0 {line, continuum) 


1215-1380 

line-by-line 

19.43 

15.89 

H 2 0 {line, continuum) 


Parameterization 

19.78 

15.67 

H 2 0 {line, continuum) 


1380-1900 

line-by-line 

30.51 

7.56 

H 2 0 (line) 


Parameterization 

30.68 

7.61 

H 2 0 (line) 


1900-3000 

line-by-line 

3.04 

5.13 

H 2 0 (line) 


Parameterization 

3.17 

5.00 

H.O (line) 


Total 

line-by-line 

Parameterization 

339.59 

342.37 

293.99 

293.03 




* The absorption due to 0 3 in Bands 4 and 6 are folded into the 0 3 absorption tables of this 
band , see Equation (4.18). 
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Table 15. Same as Table 14, except for a clear sub-arctic winter atmosphere. 


Spectral Band {cm 1 ) 

F 1 

r sfc 

/r T 

1 top 

Absorber 

0- 340 

line-by-line 

40.35 

32.21 

H 2 0 {line) 

Parameterization 

40.25 

32.08 

H 2 0 {line) 

340- 540 

line-by-line 

47.34 

52.20 

H 2 0 {line) 

Parameterization 

47.16 

51.97 

H 2 0 {line) 

540- 800 

line-by-line 

53.54 

51.13 

H 2 0 {line, continuum ), C0 2 

Parameterization 

51.73 

51.40 

H 2 0 {line, continuum), C0 2 

800- 980 

line-by-line 

1.46 

32.85 

H 2 0 {line, continuum ), 0 3 

Parameterization 

1.61 

32.84 

H 2 0 (line, continuum ) 

980-1100 

line-by-line 

3.13 

11.01 

H 2 0 (line, continuum ), O { 

Parameterization* 

3.37 

10.59 

H 2 0 (line, continuum), 0 3 

1100-1215 

line-by-line 

1.06 

9.97 

H 2 0 (line, continuum), 

Parameterization 

0.97 

10.25 

H 2 0 (line, continuum) 

1215-1380 

line-by-line 

4.61 

8.62 

H 2 0 (line, continuum) 

Parameterization 

4.86 

8.60 

H 2 0 (line, continuum) 

1380-1900 

line-by-line 

10.10 

5.01 

H 2 0 (line) 

Parameterization 

10.19 

5.09 

H 2 0 (line) 

1900-3000 

line-by-line 

0.42 

1.37 

H 2 0 (line) 

Parameterization 

0.44 

1.38 

H 2 0 (line) 

Total 

line-by-line 

162.01 

204.37 


Parameterization 

161.09 

204.21 



* The absorption due to 0 3 in Bands 4 and 6 are folded into the 0 3 absorption tables of this 
band , see Equation (4. 1 8). 
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There are a total of 13 minor absorption bands included in the LW code (see Table 1). 
These absorption bands are important in studying the greenhouse climatic effect due to human 
activities. Table 16 shows the effect of these bands on the fluxes at the top of the atmosphere and 
the surface. It is the difference between the flux with all the minor bands included and the flux 
with the absorption due to a given trace gas excluded. The concentration of these trace gases used 
is 0.2K ppmv for MO, 1.75 ppmv for CH 4 , 0.3 ppbv for CFC11, 0.5 ppbv for CFC12 , 0.2 ppbv for 
CFC22 . and 350 ppbv for C0 2 . The effect of these bands is a reduction of the outgoing longwave 
radiation al the TOA by 2.5 - 5.0 Wm 2 and an enhancement of the downward longwave radiation 
at the surface by -3 Wm 2 . 

Table 16. Effects of the minor absorption bands on the fluxes at the top of the 
atmosphere A F\ op and the surface A F l sfc as calculated using the LW radiation 
parameterization. See the text for the gas concentrations used in flux calculations. Units 
of the fluxes are Wm 2 . 



ch 4 

N : 0 

CFC's 

Minor C0 2 

Total 



Mid-Latitude Summer 



-2.22 

-1.83 

-0.47 

-0.38 

-4.91 

afV 

+0.86 

+0.58 

+0.39 

+0.86 

+2.68 



Sub- Arctic Winter 


af\ iip 

-1.00 

-1.10 

-0.19 

-0.07 

-2.36 

afV 

+ 1.21 

+1.23 

+0.39 

+0.37 

+3.20 
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11. CONCLUDING REMARKS 


This technical memorandum documents the LW radiation code developed at the Climate 
and Radiation Branch, NASA/Goddard Space Flight Center, which has been implemented in a 
number of cloud, weather, and climate models. This code includes all the major absorbers and 
most of the minor absorbers. Cloud optical properties are computed as functions of the ice/liquid 
water content, the effective particle size, and spectral bands. Aerosol optical properties are 
specified as input to the code. Scattering due to clouds and aerosols is included by scaling the 
optical thickness. To enhance the accuracy and speed of flux calculations, different approaches 
have been applied to different spectral bands for transmission calculations. In the vertical 
integration of fluxes, the adjacent layers to a level where fluxes are computed are treated 
separately in order to achieve a high degree of accuracy. Sample calculations show that this code 
can compute accurately fluxes and cooling rate within 1% and 0.15 C/day (0.25 °C/day in the 
stratosphere), respectively, of the line-by-line calculations. 

Depending upon the purpose of the LW flux calculations and the nature of the atmospheric 
model, there are a number of options available for efficiently running this code; 

1. When the upper stratosphere is not crucial in a study, one can choose the option high=. false. 
In this case, cooling rate in the upper stratosphere above the 10 ~hPa level is greatly 
underestimated. When one chooses high=.true., the cooling rate is computed accurately from 
the surface to the 0.01 -hPa level, but the computation is slower. 

2. When the absorption in the minor bands (N 2 0, CH 4 , CFOs , and minor C0 2 ) are not important 
in certain studies, one can choose trace=.false. The absorption due to those bands is not 
included. 

3. When a vegetation layer is added, the emission and reflection of both the ground and the 
canopy are included in computing the radiation emitted and reflected by the combined 
vegetation-ground surface. Set vege=.true. if a vegetation layer is included. 

4. For many cloud and meso-scale models, the size of a grid box is small. Factional cloud cover 
of a layer is either 0 or 1. Treatment of overlapping of clouds at different layers is not an 
issue, and the option overcast=.true. can be used to reduce computing time. 
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5. The cloud optical thickness either can be computed from the cloud ice/liquid water content 
and the effective mean particle size or can be an input to the LW code. Set cldwater=.true. for 
the former and cldwater=.false. for the latter, 

6. Set aerosol=.true. if aerosols are included in computing atmospheric transmission functions. 
Otherwise set aerosol=.false. 

There are still rooms for improvement to the code especially the parameterization for cloud 
optical thickness and the overlapping among cloud layers. To be available for use by atmospheric 
researchers, we have put this code and the documentation at 

http: //climate, gsfc.nasa.gov/-chou/clirad_lw. They will be updated with new developments of the 
parameterization. 
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